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ABSTRACT 

Synthetic  Aperture  Radar  (SAR)  senseis  have  many  advantages  over  electro-optical  sensors  (EO)  for  target  recognition 
applications  ,  such  as  range- independer it  resolution  and  superior  poor  weather  performance.  However,  the  relative 
unavailability  of  SAR  data  to  the  basic  research  community  has  retarded  analysis  of  the  fundamental  invariant  properties  of 
SAR  sensors  relative  to  the  extensive  invariant  literature  for  EO,  and  in  particular  photographic  sensors.  This  paper  develops 
the  basic  geometric  imaging  transformation  associated  with  SAR  from  first  principles,  and  then  gives  an  existence  proof  for 
several  geometric  scatter  configurations  which  give  rise  to  SAR  image  invariants. 
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1.  INTRODUCTION 

Current  research  into  automatic  target  -ecognition  using  SAR  has  tended  to  focus  on  the  photometric  as  opposed  to  geometric 
properties  of  SAR  sensors.  Template  based  approaches  wherein  recognition  is  obtained  by  matching  the  unknown  target 
image  chip  :o  a  precomputed  reference  image  are  popular,1’2  but  suffer  from  a  lack  of  scaleability  when  a  large  number  of 
target  classes  or  extensive  variations  in  target  appearance  are  encountered.  This  is  a  fundamental  limitation  of  such  view- 
centered  methods.3 

An  alternative  to  view  based  approaches  1o  automatic  target  recognition  (ATR)  is  an  object-centered  approach.  In  the  object- 
centered  approach,  recognition  is  accomplished  by  comparing  sensed  data  to  stored  three  dimensional  (3D)  models  that 
ideally  capture  all  of  the  essential  signature  variability  of  the  target  as  imaged  by  a  given  sensor.4  An  example  of  such  an 
approach  for  SAR  ATR  is  embodied  in  the  moving  and  stationary  target  acquisition  and  recognition  or  MSTAR  system.5  In 
MSTAR,  predicted  peak  and  region  features  are  matched  to  SAR  image  extracted  peaks6  and  regions  in  an  iterative  fashion 
to  allow  for  variation  in  the  sensed  object's  signature.  To  avoid  the  computational  cost  incurred  by  attempting  on-line 
calculation  of  radar  scattering,  MSTAR  utilizes  a  data  compression  technique7  to  extract  3D  scattering  centers  from  predicted 
xpateh-t  syrthetic  SAR  data.  Although  the  MSTAR  approach  has  been  successful  in  overcoming  some  of  the  scaleability 
limitations  of  the  template  approach,  it  does  not  exploit  all  of  the  information  inherent  in  the  explicit  consideration  of  3D 
models  and  their  projection  into  a  two  dimensional  (2D)  image. 

There  has  been  considerable  interest  in  the  object-centered  approach  to  ATR  for  optical  sensors,  and  in  particular  for 
photographic  sensors.  Numerous  examples  can  be  found  in  the  recent  image  understanding  literature.9,10,11,12  An  object- 
cen(ered  technique  of  particular  interest  is  geometric  invariance.13  Simply  stated,  a  geometric  invariant  is  an  image 
measurable  property  of  a  3D  configuration  that  is  unaffected  by  imaging  viewpoint.  Since  it’s  introduction  to  computer  vision 
in  1988, 14  geometric  invariance  has  been  extensively  developed  for  vision  applications.15  Geometric  invariants  can  be 
derived  for  single  sensor  views  of  a  3D  object  or  as  a  relationship  across  multiple  views.  One  can  term  these  two  types  of 
invariant  monoscopic  and  stereoscopic  invariants  respectively.  For  any  given  imaging  sensor,  there  are  generally  a  richer 
array  of  stereoscopic  invariants  available  than  monoscopic  invariants,  but  this  numerical  advantage  is  offset  by  the  need  for 
image  correspondence  in  order  to  use  stereoscopic  invariants.  Geometric  invariants  can  be  absolute,  or  completely 
independent  of  the  particular  imaging  geometry  (i.e.,  independent  of  sensor  position  and  orientation  relative  to  the  object),  or 
they  can  be  relative  to  a  particular  imaging  geometry.  Relative  invariants  are  typically  of  mathematical  interest  and  are  not 
generally  useful  for  object  recognition  applications,  except  to  the  extent  that  they  can  be  used  in  combination  to  generate 
absolute  invariants.  Geometric  invariants  can  be  powerful  features  for  object  recognition  since  they  alleviate  any  need  to 
search  for  a  match  over  object  pose.  To  date,  little  of  geometric  invariant  theory  has  been  applied  to  radar  frequency  (RF) 
sensors.  Payton  and  Barrett  have  examined  invariant  structures  in  stereo  HRR  sensing,16  and  Binford  has  examined  persistent 
scattering  of  BTR70  vehicle  road  wheels  in  xpatch-t  data  as  an  invariant  indexing  mechanism.17  Stuff  is  exploring  invariant 
constraints  for  application  to  3D  SAR  sensor  signal  processing,  but  little  of  his  work  has  appeared  in  the  literature  as  yet.18 

One  potential  extension  of  the  MSTAR  approach  is  to  consider  the  geometric  effects  of  the  SAR  sensor  transformation  in 
conjunction  with  the  3D  scattering  center  concept.  By  addressing  the  algebraic  nature  of  the  SAR  imaging  transformation  it 
becomes  possible  to  determine  configurations  of  3D  scattering  points  that  give  rise  to  monoscopic  SAR  geometric  invariants. 
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The  development  of  the  basic  mathematics  of  monoscopic  SAR  geometric  invariants  is  the  primary  thrust  of  this  paper. 

2.  A  HIGH  RANGE  RESOLUTION  (HRR)  RADAR  INVARIANT 

Before  developing  SAR  geometric  invariants,  it  is  useful  to  first  consider  the  one  dimensional  (ID)  radar  case.  This  allows 
for  g  simpler  analytic  development,  and  the  result  naturally  extends  to  the  2D  SAR  case.  This  also  happens  to  be  the 
historical  order  in  which  the  research  progressed.  Figure  1  gives  the  basic  imaging  geometry  for  the  HRR  sensor. 


Figure  1.  HRR  geometric  transformation 


In  this  paper,  we  shall  assume  that  any  time  mention  is  made  of  a 4  3D  point’  we  are  referring  to  a  scatterer  whose  phase 
center  is  located  at  that  point  in  3D  space,  and  not  just  any  arbitrary  3D  location  devoid  of  content.  That  is  to  say,  we  are 
focusing  our  attention  on  sources  of  RF  scattering  that  can  be  placed  at  a  particular  3D  location. 

We  approximate  the  HRR  sensor  transformation  of  a  point  in  3D  space  as  the  perpendicular  projection  of  that  point  onto  a 
vector  along  the  radar  look  direction.  Note  that  all  points  in  any  plane  perpendicular  to  the  look  direction  are  projected  into 
the  same  point  in  the  HRR  return.  This  is  an  approximation  because  the  true  transformation  projects  all  points  of  equal 
distance  from  the  radar  (within  the  radar  main  beam)  onto  the  radar  look  direction,  so  that  the  plane  in  our  approximation 
should  actually  be  a  spherical  surface.  Our  approximation  is  therefore  equivalent  to  making  a  plane  wave  assumption  for  the 
radar  wave  front  at  the  scale  of  typical  target  sizes.  In  a  typical  HRR  return,  the  signal  is  divided  into  many  discrete 
resolution  cells  that  represent  the  coherent  addition  of  all  scatters  within  a  range  bin  (in  effect  all  scatterers  in  a  3D  planar 
slab  are  projected  and  summed  coherently  into  a  single  HRR  range  cell).  We  assume  that  the  RF  scattering  is  such  that  for  at 
least  some  cells  a  single  scattering  point  in  3D  dominates  it’s  corresponding  HRR  return  cell. 


Figure  2.  A  typical  HRR  return  of  an  aircraft  object 


A  typical  HRR  return  is  shown  in  figure  2.  Since  HRR  is  a  one  dimensional  signal,  it  is  apparent  that  monoscopic  invariants 
can  only  be  obtained  as  some  function  of  the  distances  between  HRR  return  peaks.  Since  invariants  are  inherently 
dimensionless,  it  is  reasonable  to  consider  ratios  of  HRR  inter  peak  distances  as  a  likely  invariant  structure.  Our  strategy  is  to 
determine  if  there  exists  a  3D  arrangement  of  scatterers  that  would  give  rise  to  an  invariant  distance  ratio  in  the 
corresponding  ID  HRR  projection  of  the  scattering  centers.  From  figure  1,  the  projection  g  of  a  3D  point  F  given  in 

^  ^ 7  ^  coordinates  into  the  ID  HRR  signature  is  easily  found  given  the  radar  look  direction  unit  vector  ^by 

where  o  is  the  arbitrary  origin  of  the  HRR  coordinate  system.  Note  that  without  loss  of  generality  a  can  be  assumed  to 
include  any  scale  factors  needed  for  the  HRR  down  range  variable.  This  gives  the  fundamental  HRR  transformation 
projection  constraint 

a-ffj,(r)  =  {r-o)-s  n , 
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We  are  interested  in  determining  conditions  under  which  we  have  distance  ratio  invariants  in  HRR.  This  generally  implies  at 
least  three  distinct  points  in  3D  space.  We  denote  distinct  3D  points  by  P'  ,  3D  Euclidean  distances  between  any  pair  of  two 

distinct  points  ^  by  ^  ^ ,  and  the  corresponding  ID  HRR  distances  by  K  a*  I .  We  can  then  relate 

the  3D  and  ID  distances  between  distinct  points  by 

a,  "a,  =(fe  -®))  >=  (* 


where  Pfl  is  the  unit  direction  vector  from  to  P* .  Therefore 

Note  that  an  added  benefit  of  considering  inter  peak  distances  in  the  HRR  signature  is  the  elimination  of  the  need  to  consider 
the  HRR  coordinate  origin,  so  that  we  immediately  gain  invariance  to  3D  translation.  A  ratio  of  distances  between  dominant 
scatterers  on  the  HRR  line  is  related  to  the  corresponding  3D  distances  between  the  scatterers  by 

4  (a*  >aJ  d3 (p* ,pt )  |pM  •  (5) 


Therefore  HRR  distance  ratios  give  an  invariant  if 


| PjlA  = 

\p*/\ 


constant 


regardless  of  the  orientation  ^of  the  radar  look  direction.  Clearly  a  sufficient  condition  for  this  to  arise  is  for  P*'  and  Pfl  to 

ri  —  +  rj 

be  parallel  vectors,  that  is,  ^}t  *9L ,  and  in  this  case  the  constant  is  unity  and  we  have  an  absolute  invariant.  For  arbitrary 

^it  is  readily  seen  that  this  must  also  be  a  necessary  condition  as  well.  This  can  be  shown  as  follows.  For  ease  of  notation 
denote  ^  and  Pf*  by  P  and  ^  respectively.  Our  constraint  then  is 

tp * M  .i ~  ~ 

pr^r  =  f\P*q*A  =  constant 

^  ,(7) 

and  we  are  interested  in  the  relationships  between  P  and  ^  satisfying  this  constraint  for  arbitrary  subject  to  and 

so  that  this  constraint  does  not  become  undefined  or  trivially  satisfied  respectively.  Let  ^  ^  .  Then  since  ^  ,  we 

have  ~  * .  Similarly  for  ^  P  we  obtain  fiP’^’P)  ~  * .  Therefore  for  arbitrary  ^  ^  ^ s  case 

]p  -[$'$  =0  =  ±r>(g) 

Since  this  equation  must  hold  for  arbitrary  *  ,  it  holds  in  particular  for  ^  ^  and  so  we  obtain  P  ^  ~ 1  or  P  “  ,  i.e., 

parallel  vectors.  This  proves  that  parallel  vectors  is  both  a  necessary  and  sufficient  condition  for  an  absolute  ratio  distance 
invariant  in  HRR.  Note  that  Payton  and  Barrett16  derive  a  similar  result  for  the  case  of  three  collinear  points  as  a  special  case 
of  a  two  view  stereoscopic  HRR  invariant  for  four  coplanar  scattering  centers,  which  they  term  a  "Velten  Invariant."  As  the 
name  might  imply,  this  invariant  was  first  suggested  by  the  author,  and  Payton  and  Barrett  examined  HRR  invariants  at  the 
author’s  request  under  contract  to  AFRL. 

3.  SAR  GEOMETRIC  INVARIANTS 


We  begin  as  we  did  for  the  HRR  case  by  deriving  an  approximation  for  the  SAR  3D  to  2D  imaging  transformation.  As  in  the 
HRR  case,  this  is  an  approximation  of  the  SAR  imaging  projection,  as  the  true  physical  situation  the  projection  is  also  along  a 
spherical  front  centered  at  the  synthetic  aperture  phase  center.  However,  for  objects  of  limited  extent  such  as  vehicles,  the 
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error  introduced  by  assuming  a  perpendicular  projection  into  the  slant  plane  is  negligible,19  and  we  are  in  effect  assuming 
that  the  point  is  in  the  far  field  of  the  synthetic  aperture. 


Figure  3.  SAR  imaging  geometry 


Denote  3D  world  coordinates  as  ix^,z) ,  general  (column)  vectors  in  3D  world  coordinates  by  f  ,  and  unit  vectors  in  world 
coordinates  by  r  .  where  r  is  any  Roman  symbol.  The  SAR  slant  plane  position  and  orientation  are  given  by  unit  normal 
vector  n  located  at  ^  in  the  3D  world  frame.  We  specify  a  2D  coordinate  system  lying  in  the  slant  plane  (i.e.,  the  SAR  image) 

by  perpendicular  vectors  u  and  v  5*v  ,  where  s*  and  are  slant  plane  coordinate  scale  factors  and  **  and  ^  are  unit 
vectors.  Finally,  we  denote  locations  in  the  2D  SAR  image  coordinate  system  using  Greek  symbols  e.g.  a  ,  where  a  is  an 


ordered  pair  specifying  components  in  the  '  ’ 
generally  the  result  of  a  dominant  3D  scatterer. 


U ' directions.  As  in  the  HRR  case,  we  will  assume  that  SAR  image  peaks  are 


Approximate  the  SAR  image  of  a  point  at  P  as  the  perpendicular  projection  of  P  into  the  slant  plane.  We  wish  to  determine 
the  SAR  image  coordinates  of  P ,  or  a  in  figure  3,  as  a  function  of  the  know  slant  plane  geometry.  Clearly  a  is  the  slant 
plane  perpendicular  image  of  P  %  .  Since  a  ~  P  ^  is  a  vector  whose  origin  is  at  the  slant  plane  origin,  we  can  proceed  by 

finding  the  change  of  coordinate  matrix  T  than  transforms  vectors  in  coordinates  with  their  origin  at  %  into 


(u,v,n) 


^  7  *  coordinates.  Clearly  1 

can  then  easily  determine  a  from 


T  =  [u  v  w] r 


,  so  that  a  transforms  to 


\'=Ta  =  {u-a,va,na)r  (w) 


system.  We 


&=  Pa1 


1  0  0 
0  1  0 


1  0  0 


o  i  or*  o  i  or* 


1  0  0 


ur  um  u 


u  u  it 


*  '  la  -  (u  a,v-  a) 

V.  V.  V  '  ’  / 


iti  IT 

Note  that  we  can  express  T  as  the  product  of  scaling  and  rotation  matrices  1 1  and  imf  respectively. 


s  0  0  QT 

U.  || 

T=TJ^=  0  *.  0 

0  0  1  nT 


Therefore  the  SAR  image  coordinates  a  of  a  3D  point  at  P  are 


=[» 
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For  the  case  of  a  SAR  system  with  ‘square’  pixels,  i.e.,  for  which  the  along  and  cross  range  image  dimensions  are  of  the  same 
resolution,  we  have  u  v  and  the  SAR  image  coordinate  expression  becomes 

^  i?rr  "^)s  sM(^  *) 


We  shall  only  consider  SAR  systems  with  square  pixels  in  this  paper.  It  is  instructive  to  consider  the  overall  algebraic 
structure  of  the  SAR  imaging  transformation.  Prior  to  projection  into  the  slant  plane  we  have 

This  transformation  is  a  special  case  of  an  affine  3D  to  3D  transformation,  since  the  matrix  is  taken  from  the  special 
orthogonal  group  (i.e.,  othronormal  matrices)  of  dimension  3  over  ^ ,  which  is  a  subgroup  of  the  general  linear  group 

^  invertible  matrices  over  .  This  transformation  is  a  group  action  in  3D  space.  To  get  the  SAR 
transformation  this  group  action  is  followed  by  a  direct  subspace  projection  using  the  projection  matrix 

'=[■  °  i 

t>  l  o 


into  the  slant  plane  (or  ^  )  as  given  above.  Note  that  by  use  of  homogenous  coordinates20  we  can  express  the  affine 

=  T 

transformation  of  a  point  ”  via  a  single  homogenous  4x4  transformation  matrix  A  as 


=  A  -r 


v  -  v  q  p 

«r  ~nq  1 

0  0  0  1 


and  the  SAR  transformation  can  be  expressed  as  a  2x4  matrix  operating  on  3D  points  expressed  in  homogenous  coordinates 
via 

/:»*-»  »*:/(?)  =  5,  Z?=s(“r 


Thus  there  are  6  degrees  of  freedom  for  the  SAR  transformation  needed  to  specify  the  location  (3  degrees  of  freedom  in  ®  ), 

7* 

orientation  (2  degrees  of  freedom  in  the  rotation  matrix  ™r ),  and  the  coordinate  scaling  (1  degree  of  freedom)  in  the  slant 
plane.  Note  that  the  translation  component  of  the  SAR  transformation  matrix  will  cancel  out  for  line  segments 

/  v‘  -V \  (?]"[?!  =  s[^ 


making  it  convenient  to  work  with  distances  between  points  instead  of  individual  points  since  by  doing  so  we  immediately 
obtain  3D  translation  invariance.  Now  that  we  have  derived  the  essential  transformation  information,  we  can  explore 
invariant  constraints.  We  can  easily  relate  the  distance  between  any  two  distinct  points  in  the  SAR  image  to  their  3D  distance 
as 

d*  ( ^  ^ ) s  IN  ■  ^  II = IM*  _a)I= Ik^  (pi  -  p*  )ai  II = <*>  b.  >  Pt  JlkA.  II  ne. 


where  P*1  is  the  unit  vector  in  the  P'  P’1  direction,  and  ^  is  Euclidean  i 
distances  between  2  distinct  pairs  of  SAR  image  points 

,  j| )  _  dj  (ft.ft llkaill 

d,  (ft  .A  1^1 


is  Euclidean  distance  in  k  dimensions.  Therefore  a  ratio  of 
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will  be  invariant  to  the  SAR:  imaging  geometry  provided 

f  *  .  fMpoJ 

1=  constant 

II  Pvt  .{20) 

Since  we  seek  3D  point  configurations  that  will  satisfy  this  constraint  for  arbitrary  it  is  obvious  that  one  such  condition  is 

that  the  distinct  point  pairs  form  parallel  lines  in  3  space,  i.e.,  “  -Pii .  Under. this  condition  we  would  have  an- absolute 

invariant  and  the  constant  would  be  unity  in  this  case.  Whether  this  sufficient  condition  for  an  absolute  invariant  is  also  a 
necessary  condition  remains  to  be  determined. 


We  shall  show  that  (20)  is  also  a  necessary  condition  for  absolute  invariance.  To  prove  this,  we  shall  proceed  as  above  for  the 
HRR  case.  For  ease  of  notation,  replace  ^  and  with  "P~ and  ^respectively.  Then  we  have 


[tL-J?  Hv  pf 


and  we  are-  interested  in  bounds  on  ^  for  arbitrary  ^ancf  ?  .Note  that  botLlhc  numerator  and  denominator  of  this  expression 

lie  between  0  and  1,  since  each  is  the  magnitude  of  the  projection  of  a  3D  unit  vector  into  a  2D  subspace.  Now  let  u  ~  P . 
Then 

/  -  -  ;  ■  ‘  it 

MWW4  -rrr. 


For  U  ^  we  obtain 


Therefore  for  arbitrary  “  and  $  t  f  *  5  which  na^urn,  implies 


Let  U~F.  Then 


Let  u  ~  % .  Then 


(£  •  pf  +  {v  ■  pf  =  (u  qf  +  (v -qf 


*{v  pf  ‘=  i 


Hence  for  arbitrary  u 


(v  •  qf  -  (v  •  pf  =»  (v  ■  q  -  v  •  ,p)(v  q  +  v-jfy  =  0=*v'q  =  ±vjp 


A  __  —  +yi 

and  since  v  is  also  arbitrary,  we  obtain  “  ~  or  parallel  vectors.  Hence  parallel  vectors  is  a  necessary  and  sufficient 
condition  for  an  absolute  ratio  distance  invariant  in  SAR  imagery. 


3.1  Generalizations  of  the  parallel  vector  invariant 
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Having  established  the  invariance  of  parallel  vectors  under  HRR  and  SAR  imaging  geometry,  a  generalization  of  this  result 
immediately  suggests  itself.  It  is  clear  that  two  triangular  areas  (i.e.,  two  triples  of  co-planar  scattering  centers)  lying  in 
parallel  planes  also- project  invariantly  (figure  4). 


Figure  4.  SAR  projection  of  triangular  areas  lying  in  parallel  planes 


This  is  because  the  area  of  a  triangle  is  proportional  to  the  magnitude  of  the  cross  product  of  any  two  of  it’s  sides  expressed 
as  vectors  in  3  space  (figure  5).  Since  the  cross  product  is  perpendicular  to  it’s  base  vectors,  then  any  cross  product  pairs 
taken  from  two  triangular  areas  lying  in  parallel  planes  will  also  be  parallel,  and  by  the  results  of  the  previous  sections,  the 

be  preserved  under  either  the  SAR  or  HRR  transformation. 

j l*|k  =  j|a(b|sin^-||ff  xi| 


n 

Figure  5.  Area  of  triangle  formed  by  vectors  &  and  & 

Since  the  SAR  transformation  substantially  preserves  the  scattering  center  geometry,  this  observation  may  prove  useful  as  a 
feature  in  SAR  imagery.  Thus  there  would  be  a  rich  set  of  potential  invariants  in  the  case  of  two  parallel  planes  of  coplanar 
scattering  centers,  a  situation  that  may  be  common  in  SAR  imagery  of  vehicles  and  other  man-made  structures.  Note  that  4 
coplanar  points  with  no  three  colinear  would  also  give  rise  to  invariants  as  a  special  case  of  the  above,  since  there  would  be  4 
distinct  triangles  of  scattering  centers.  From  4  distinct  triangles  6  independent  ratio  invariants  can  be  formed. 


2.  SAR  quasi-invariants 


In  many  cases  it  may  not  be  wise  necessary  or  even  to  demand  an  absolute  invariant  to  obtain  a  feature  useful  for  recognition 
in  SAR  imagery.  This  is  because  the  fluctuation  in  scattering  center  amplitude  and  phase  center  location  as  a  function  of  SAR 
imaging  geometry  is  such  that  scattering  peak  location  persistence  is  limited  to  a  few  degrees  in  most  cases.21  It  is  therefore 
prudent  to  examine  whether  a  more  general  geometric  configuration  can  be  found  that  provides  feature  stability  with  pose  in 
a  local  sense.  Such  local  invariants  are  termed  quasi-invariants  by  Binford.22 


Following  Binford’s  definition  of  quasi-invariant,  to  determine  the  conditions  under  which  ratios  of  distances  in  a  SAR  image 
give  rise  to  a  quasi-invariant,  we  expand  our  expression  for  a  SAR  image  distance  ratio  as  a  multivariate  Taylor  series  in  the 

imaging  transformation  variables  ^  and  ^ .  We  then  examine  conditions  on  ?  and  ^  under  which  the  coefficients  of  the  linear 

terms  vanish  in  order  to  determine  quasi-invariant  ratios.  For  the  purposes  of  this  analysis,  it  is  sufficient  to  consider  ^ 

constant  7  implies  constant  7  .  We  select  fixed  ^  and  ^  at  ^  and  ^respectively.  Thus 

/i  =  =/IU+v^,‘fiLils‘^l+v^I^LI?"i?*l+0  ^ if *)  (2g) 


since 
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We  are  interested  in  determining  ?  and  ^  that  satisfy 


that  4s,  finding  P  and  ^  that  force  the  coefficient  of  the  linear  term  in  the  Taylor  expansion  to  vanish.  For  ease  of  notation  we 
drop  the  0  subscript  on  u  and  v  .  Note  that  since  ^  is  symmetric  with  respect  to  “  and  ^  ,  then  ^  will  also  be  symmetric 
with  respect  to  ^  and  v  .  Thus  we  will  get  the  same  constraint  result  on  P  and  ^  from  forcing  either  linear  coefficient  to 


vanish.  Hence 


vj2£  = 


♦(*•£)  pf_  +(*  tO(u' jte 

[M1*  (*•$)*) 

^  (fcg)2  + ~  (fc  •  £)*  fog/ 

(M*  +  (v  ^)s ) 

„  fasffe  pY  -(»■ -p)*  (“•«)* 


( 

[(*•«) 

(“«)'  *1 
M)'  - 

r  *«)* ) 
(E-p)1 

l 

L 

[«•*)/ 

i 

;mm 

2 

Note  that  the  denominator  cannot  vanish  for  non-trivial  ^  provided  ^  ^  .  In  a  similar  manner  the  denominator  of 

^  will  not  vanish  provided  P  *  .  Note  that  the  case  of  P  and  ^  both  parallel  to  ^  would  obviously  give  an  invariant 

for  any  choice  of  n  ,  so  we  can  safely  assume  a  non-trivial  denominator  for  at  least  one  of  the  linear  terms  in  the  Taylor  series 
expansion.  Equation  (30)  then  reduces  to 

[vqup-v-puq)(vqu-p+v‘puq)=d 


and  thus 


vqup  -  tvpirq 


If  we  assume  that  none  of  the  terms  vanish  then 


UL  =  ±v'p 

“«\P  .(33) 


Note  that  each  side  of  this  equation  can  be  interpreted  as  the  slope  of  the  projection  of  a  three  space  vector  (P  or  9  )  into  a 

2D  subspace  (the  ^  ?  *  plane).  Hence  the  geometric  interpretation  of  this  condition  is  that  the  projections  of  ^and  ^  into  the 
slant  plane  are  either  parallel  (+),  or  that  they  possess  mirror  symmetry  in  slope  (-),  e.g.,  that  have  slopes  of  equal  magnitude 
but  opposite  sign.  We  can  therefore  state  that  point  sets  that  generate  line  segment  pairs  with  slopes  of  equal  magnitude  when 

projected  into  the  slant  plane  are  a  sufficient  condition  for  a  geometric  quasi-invariant.  Note  that  the  ®  component  of  P  and 

^  are  unconstrained  in  this  case.  The  strength  of  this  quasi-invariant  (that  is,  the  viewing  region  over  which  it  remains 

approximately  constant)  would  clearly  depend  on  the  relative  magnitudes  of  the  ”  components  of  P  and  ^ ,  as  a  large 
B  component  would  give  the  feature  more  ‘wobble’  as  the  imaging  geometry  is  varied  in  azimuth  about  the  target. 

4.  CONCLUSIONS  AND  FUTURE  WORK 


http://www.mbvlab.wpafb.af.mft/papers/6tiws/geohrr.htrnr 


4/29/04 


Distance  Ratio  Invariants  in  SAR  Imagery  - 


Page  9  of  1 1 


% 

This  paper  introduced  the  concept  of  pose  invariant  features  for  object  recognition  in  SAR  imagery,  and  showed  some  basic 
3D  geometric  arrangements  that  give  rise  to  such  features.  The  results  show  that  at  least  theoretically,  parallelism  and 
coplanarity  of  scattering  center  locations  gives  rise  to  potential  SAR  image  geometric  (absolute)  invariants.  It  remains  to  be 
seep  whether  these  configurations  provide  a  sufficTeirtf3rrich  set  ofgeometric  relationships  to  allow  for  a  useful  number  of 
SAR  image  invariants  for  vehicle  object  recognition.  This  will  be  the  subject  of  future  research.  A  potentially  fruitful 
direction  is  to  consider  the  use  of  general  invariant  constraints  thatarise  between  affine  3D  to  3D  transformations  and  affine 
2D  to  2D  transformations  for  3D  point  sets  and  their  corresponding-2D  point  set  images  under  the  action  of  the  SAR  image 
projection  operation.  This  would  be  similar  to  researchr  being  done  by Weiss  for  the  more  general  case  of  the  affine  and 
projective  optical  camera.23,24,25 

The  MSTAR  program  has  recently  released  a  simplified  version  of  a  3D  scattering  center  based  SAR  synthetic  target 
signature  image  prediction  tool  for  T72,  BMP2,  and  BTR70  vehicles  called  Predict  Lite.  Using  the  Predict  Lite  3D 
scattering  center  signature  representation  of  these  vehicles,  it  will  be  possible  to  analyze  3D  scattering  center  configurations 
of  a  given  target  object  to  determine  potential  SAR  image  geometric  invariants.  The  utility  of  such  potential  invariants  can  be 
assessed  by  examining  2D  peaks  in  Predict  Lite  synthetic  imagery  generated  from  the  corresponding  3D  scattering  center 
data  base.  If  the  synthetic  imagery  suggests  that  the  invariants  have  recognition  discrimination  power,  then  the  invariants  can 
be  evaluated  using  measured  SAR  imagery  of  tactical  vehicles  recently  released  by  D  ARP  A/ISO  and  AFRL. 

As  part  of  the  research  plan  outlined  above,  it  will  also  be  necessary  to  assess  the  impact  of  certain  signal  processing  effects 
in  the  SAR  image  formation  process  on  invariant  stability  and  accuracy.  Major  SAR  image  formation  effects  likely  to  be 
important  are  the  radar  2D  impulse  response,  scintillation,  discretization,  and  antenna  weighting.  These  effects  are  similar  to 
impulse  response  and  radial  distortion  effects  in  optical  sensors.  Experience  has  shown  that  accurate  calculation  of  optical 
image  invariants  requires  compensation  for  such  optics  distortions,  so  it  is  reasonable  to  expect  that  the  analogous  situation 
will  arise  in  SAR  imagery.26 

Since  the  SAR  sensor  has  a  finite  2D  impulse  response,  scattering  energy  that  would  be  concentrated  at  a  point  in  an  ideal 
infinite  resolution  sensor  is  instead  spread  out  over  the  resulting  image  in  the  well  known  sin(x)/x  pattern.  This  has  obvious 
practical  implications  for  the  stability  of  image  invar iants  derived  from  the  3D  scattering  center  representation  in  Predict  Lite. 
Fortunately  Predict  Lite  includes  the  SAR  impulse  response  as  part  of  the  synthetic  image  generation  process,  so  that  the 
impact  of  this  phenomenon  can  be  assessed  using-synthetic  data;  Antenna:  weighting  functions  have  an  impact  similar  to  the 
effects  of  the  2D  impulse  response. 

Scintillation  is  a  phenomena  that  arises  as  a  consequence  of-the  finite  spatial  resolution  in  SAR.  A  given  range  and  azimuth 
bin  may  contain  contributions  from  many  small,  unresolved  scatters  which  when  coherently  summed  together  cause  a  net 
response  that  is  highly  variable  with  small  changesTn  the-SAR  imaging-geometry  (i.e.,  the  6  parameters  of  the  SAR  imaging 
transformation  from  section  3).  This  will  tend  to  make  some  3D  scattering  centers  project  in  an  unreliable  manner  into  the 
SAR  image. 

The  fact  that  a  SAR  image  is  a  discrete  representation  of  an  undcrly  ing  continuous  function  will  effect  the  accuracy  and 
stability  of  invariant  values.  Each  image  pixel  represents  a  coherent  average  of  the  underlying  scattering  energy  with  it,  and 
this  will  clearly  tend  to  shift  scattering  positions  awayfromtheiridealized2D  locations  as  predicted  by  the  basic  SAR 
projective  geometry  model.  Use  of  state-of-the-art  SAR  sub-pixel-peak  detection  algorithms  may  mitigate  these  effects 
somewhat,  but  they  cannot  eliminate  them  altogether. 
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